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ABSTRACT 

The mass of the central black hole in Seyfert galaxies and QSOs can be determined from the broad 
emission lines and the reverberation method. Using the measured black hole mass and the bolometric 
or ionizing luminosity, the accretion rate can be estimated. Compiling a sample of Seyfert 1 galaxies 
and QSOs with reliable central masses, estimated ionizing luminosities and X-ray spectral slopes, we 
find that the X-ray spectral slope strongly correlates with the accretion rate. The objects in the sample 
are found to be distributed in two distinct classes in the spectral index versus Log(Lion/iEdd) plane. 
We argue that these two classes may correspond to ADAF and thin disk accretion. The observations 
of a "two-state" Seyfert 1 galaxy, 1H0419-577, confirm our results. Detailed fitting of the spectra of 
individual Seyfert 1 galaxies and QSOs using ADAF and/or thin disk models should further clarify the 
two-class classification. 

Subject headings: accretion, accretion disks-galaxies: active-galaxies: quasars 



1. INTRODUCTION 



son & Wandel 1999). With known MBH masses, the ac- 



The QSO and active galaxy phenomena can be success- 
fuUy explained by the aGc rction of gag onto central mo s - 



cretion rates of active galaxies and QSOs can be estimated 



from their bolometric o r ionizing luminosity (Wandcl, Pc- 



bive black holes (MBHs) (Lyndon BoU 1969 h Roos 198d ) 



tcrson & Malkan 1999, hereafter WPM). Furthermore, a 



Recent progress in both theory and observation strongly 
, supports this explanation. First, compact massive dark 
objects — possibly massive black holes — have been detected 
\ in the cores of some nearby galaxies based o n high spatial 
. resolution observations of stellar dynam ics (Kormendy & 
' Richstone 1995; Magorrian et al. 1998). Second, the ra- 
diation spectra from the nuclei of Sgr A* and NGC 4258 
can be modelled well by the advection-dominated accre- 



tion fiows (ADAF) onto MBH s ( [Narayan, Yi fc Mahade 
van 1995|; |Lasota ct al. 1996| ). In addition, the iron Ka 



line profile suggests the existence of a MBH and cold disk 



mate rial in the center of Seyfert 1 galaxies (Iwasawa et al 



1996). 



Currently, four sol utions of the accreti on process around 
MBHs are known ( phen ct al. 1995|). The most fa- 
mous two are the thin disk n ioder^ Shakura fc Sunyacv 
' I973p and the ADAF model ( flchimaru 1977|; |Rees ct al, 
19821 Narayan & Yi 1994, 1995a, 1995b; lAbramowicz ei 



iifferl 



1995 ), which generally result from dift'erent accretion 
rate. The radiati on spectra output of these models are dis - 
tinctly different ( Narayan, Mahadevan fc Quataert 19981 ). 
It is therefore reasonably easy to tell from observations 
whether a system has a thin or ADAF disk. The ADAF 
model provides a good fit to the spectra of Galactic Black 
Hole X-ray Binary sources (BHXBs), which should be a 
l ow mass accretion version of Seyfert 1 galaxies and QSO s 



dEsin, McChntock fc Narayan 199^ ; |Esin et al. 1998] ) 



?his implies that at least some Seyfert galaxies and/or 
QSOs, which were believed to have a thin disk accretion, 
may instead accrete material via the ADAF mode. Thus 
the question of whether a QSO or an active galaxy has an 
ADAF or a thin accretion disk is of great current interest. 
To date, the most reliable estimate of masses of central 



MBHs in Seyfert 1 galaxies uses reverberation data (Peter 

Tech 



remarkable correlation between the soft X-ray spectral in- 
dices and the full width at half maximum (FWHM) of 
H/3 has been found by Boiler, Brandt & Fink (1996) and 
Wang, Brinkmann & Fink (1996). A possible interpreta- 
tion is that the correlation is cau sed by the variation of 



accretion rate in diff erent objects (Wandel & Yahil 1985 



Pounds et al. 1994). The change in soft X-ray spectral 
slope from hat to steep may reflect the change in accre- 
tion mode from a low state to a high state. It is reason- 
able to constrain theoretical models by investigating the 
statistical relation between the soft X-ray spectra and the 
estimated accretion rates for a sample of objects. 

In this letter, we show that a sample of QSOs and 
Seyfert 1 galaxies with available estimated MBH masses 
can be classified into two distinct classes with different 
accretion rates. These two classes should correspond to 
different accretion disks, i.e. ADAF and thin disk. In 
each class, the soft X-ray spectral slope strongly corre- 
lates with the estimated accretion rate. We describe the 
available data in §2. The statistical analysis is presented 
in §3, and an explanation of the correlation between the 
X-ray slope and the accretion rate is discussed in §4. 

2. THE DATA 

Assuming that the line-emitting matter is gravitation- 
ally bound and hence has a near-Keplerian velocity dis- 
persion (indicated by the line width) , WPM estimated the 
mass of the central black holes for 17 Seyfert 1 galax- 
ies and 2 QSOs using reverberation data. The central 
MBH masses were also measured by the same method 
for three additional Seyfert 1 galaxies (Mrk279, NGC3516 
and NGC4593) which were not included in WPM's sam- 
ple ( |Ho 199^ ). Laor (1998) estimated the MBH masses 
for 19 QSOs using the H/3 width and the empirical rela- 
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tion r |BLR.(H/3) ^ 15L'^'/ (|Kaspi et al. 199^ where vblr 
is the size of H/? emitting region in the broad line region 
and L44 = L(0.1 — 1/im) in units of 10'*''ergs s^^. Laor's 
method may overestimate the black hole mass, as WPM 
pointed out that the slope of the empirical relation may 
be fla tter than 1/2. Fortunately, there is a common ob- 



ject t( |) Laor and WPM's sample, the QSO PG0953+414, 
whose mass was estimated to be 3 x lO^M© by Laor, and 
1.5 X lO^MfT^ by WPM. We will thus use a calibration fac- 



tor of 1 0.5 for the central black hole masses of the QSOs in 
Laor's sample. 

Generally, we can estimate the accretion rates for QSOs 
and Seyfert 1 galaxies from their bolometric luminosity. 
In WPM's sample, they list the ionizing luminosity of 
indiyidual objects, derived from the reverberation and 
photoionization method. Following that method, we es- 
timate the ionizing luminosities for Mrk279, NGC3516 
and NGC4593. Laor (1998) estimated the bolometric lu- 
minosities for the QSOs in his sam ple using the empiri- 
cal relation Lboi = 8.3i/L^(3000A) ( |Laor fc Draine 1993[ ) 
for the objects of the sample. For the common object, 
PG0953-I-414, we compare the estimated Lboi in Laor's 
table 1 with Lion in WPM, and find that Lboi = S.BLjon 
and iion/'t'Edd — 0.53Lboi/iEdd- We take this relation as 
an estimation of the Eddington ratio of ionizing luminos- 
ity for the objects in Laor's sample. We list in table 1 
the central black hole masses and the Eddington ratios of 
ionizing luminosity for all the objects studied. Since the 
ionizing and bolometric luminosities have been estimated, 
and the black hole mass determination is the most reliable 
one, these values give a direct measurement of the Edding- 
ton ratio, L/L^dd cx; iboiZ-^BH, and hence the accretion 
rate m — M / M-^dd- 

In general, the X-ray emission spectra of QSOs and 
Seyfert 1 galaxies can be well fitted by a power-law {fn oc 
where a is spectral index and F = 1 + a is photon 
index). The soft X-ray spectra for most objects in table 1 
are available in the literature. We list the spectral indices 
in table 1 both for the soft X-ray range (0.1-2.4keV) ob- 
served by ROSAT and for the hard X-ray range (3-lOkeV) 
observed by ASCA. For the soft X-ray spectra, we adopt 
the spectral indices of the fitting by using a power-law 
with the column density absorption as a free parameter. 
For several objects without such fitting, we use the spec- 
tral indices determined assuming a fixed column density 
absorption as given by the Galactic column density. There 
is little difference in the ASCA spectra fitted by different 
Fe Ka line models. We adopt the photon indices given by 
Nandra et al. (1997) in their table 4 by fitting 3-lOkeV 
spectra with a power-law plus a Schwarzschild geometry 
disk line. The ASCA photon indices of several objects, 
which are not included in Nandra et al. (1997), are also 
listed in table 1. 

3. STATISTICAL ANALYSIS 

Figure 1 shows the spectral indices as a function of the 
ratio Log(iioii/-/jEdd), and hence the accretion rate. We 
notice that several objects in our sample, such as 3C273, 
3C120, 3C390.3, 3C323.1, PKS1302-102 and PKS2135- 
147, are radio- loud galaxies or QSOs. The X-ray spectra 
of radio-loud galaxies and QSOs are generaUy flatter than 
those of radio-quiet ones (Laor et al. 1997; Wang, Lu & 



Zhou 1998). It may be caused by the X-ray emission from 
the relativistic jet. Indeed, all the radio-loud objects ex- 
hibit flat soft X-ray spectra in our sample. They will be 
excluded in the following statistical analysis. In addition, 
the soft X-ray spectrum of NGC3227 is ver y flat, which 
may be ca used by a dusty warm absorber (Komassa & 



Fink 1998). It will also be excluded in the statistical anal- 



ysis. 

We apply the Pearson linear correlation test (Press 
et al. 19*9^ ) to the aRosAT versus Log(Lion/iEdd) re- 
lation for radio-quiet objects. A significant correlation 
is found, with a linear correlation coefficient R — 0.80 
(P). — 1.4 X 10^*"). The soft X-ray spectra become steeper 
when the Eddington ratio of the ionizing luminosity in- 
creases. It is evident that there is a discontinuity near 
the critical point of Log(Lion/iEdd) ~ —1.4 and the data 
in figure |l| are distributed in two separate classes. One 
is the lower Log(Lion/iEdd), representing a lower accre- 
tion rate class, the other is the higher Log(Lion/iEdd), 
representing higher accretion rate class. The two dis- 
tinct classes may correspond to the accretion of ADAF 
and thin disk since the break point corresponds to the 
critical threshold between the two models (see following 
section). The Pearson linear correlation coefficients for 
the two different classes are R = 0.88 (Pr = 1.4 x 10""*) 
and R = 0.89 [P, = 2.0 x 10"*^), respectively. They can 
be fitted by straight hues: cirosat = (2.69 ± 0.24) -|- 
(0.69 ± 0.12)Log(Lion/iEdd) for the ADAF class, and 
"ROSAT = (2.00 ± 0.06) + (0.60 ± 0.08)Log(Lion/iEdd) for 
the thin disk class. 

Most Seyfert 1 galaxies in the WPM sample belong to 
the ADAF class and all QSOs in Laor's sample belong 
to the thin disk class. It is interesting to note that five 
Seyfert 1 galaxies (MrkllO, Mrk335, Mrk509, Mrk590 and 
NGC4051) join in the thin disk class. Those Seyfert Is 
have small central black holes but normal Seyfert lumi- 
nosities. Those support the existence of the two-class dis- 
tribution. 

The relation between Fasca and Log(Lion/-^Edd) is pre- 
sented in figure |[ There are only 16 objects in figure || 
because the others have no ASCA data. We find that 
there is a linear correlation between Fasca and iion/-t'Edd 
with a correlation coefficient R = 0.74 (P^ = 0.006) (after 
we exclude the radio-loud objects and NGC3227). The 
hard X-ray spectra become steeper when Log(Lion/iEdd) 
increases. As figure || shows, the photon indice Fasca in- 
creases with increasing accretion rate for ADAF class. It is 
hard to determine the behavior of the hard X-ray spectra 
vs. Log(Lion/iEdd) for the thin disk class because there 
are only 4 objects. Further observations are needed to 
confirm the classification. 

4. DISCUSSION 

We find two distinct classes in the aRosAT versus 
Log(Lion/iEdd) plane for a sample of Seyfert 1 galaxies 
and QSOs. This classification is striking because it re- 
veals the existence of two accretion modes in these ob- 
jects. The critical threshold between the two classes is 
Log(Lion/iEdd) — —1.4 (see in figures |l| and figures ||). 
The corresponding Log(Lboi/iEdd) is about -1.1, which is 
just the critical accretion rate Log(mcrit) ^ — 1.1 (see the 
figure 7 in N arayan, Mahadevan fc Quataert 1998 ). Since 
Wcrit ~ (Esin, McClintock fc Narayan 1997), where 
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a is the viscosity parameter, the observed discontinuity in disk is Compton-scattered by the hot gas in the ADAF 



figure s [l| suggests a ~ 0.3. Below mr,,)- , the accretion is via It is similar to the intermed iate state in BHXBs (Esin, 



an ADAF at small radii and a thin disk at large radii. Be- 



McClintock & Narayan 1997). The reflection component 



cause much of the viscously generated energy is advected 
to the black hole and £hn i/-^Edd — 30m^ (for the viscosit y 



parameter a — 0.3, Esin, McClintock & Narayan 1997), 
we can derive that m of the ADAF class in our sample 
is distributed in the range (0.01-0.08), which corresponds 
to the low and intermediate states of accretion in BHXBs. 
The existence of the ADAF class in Seyfert Is and QSOs 
is also supported by the fact that NGC 4151 has almost 
an identical spectrum to the BHXB GX339-4 ( ^dziarski et] 
al. 19 j^). Above rhcriti the accretion is via a thin disk and 



the accretion rate m is approximately -/jboi/^Edd- It may 
correspond to the high and/or very h igh states in BHXBs 
(Esin, McChntock & Narayan 1997). According to the 
best fitting straight lines for the two classes, we derive that 
aRosAT on m^-^^ for the ADAF class, and OfRosAT « to^-^" 
for the thin disk class. 

For objects in the thin disk class, it is easy to understand 
that the soft X-ray slope steepens with the increase in ac- 
cretion rate and the soft X-ray spectra are steeper than the 
hard X-ray spectra. The temperature of the inner thin disk 
goes up while the accretion rate increases, therefore the big 
blue bump moves into EUV and/or soft X-ray range and 
the soft X-ray spectrum becomes steeper. As to increases, 
the disk flux to irradiate the corona increases. This will 
cause the corona to cool more efficiently because of Comp- 
ton cooling, and so the hard X-ray slope will steepen with 
the increasing to. This cannot be confirmed in our sample 
because of the lack of ASCA observations for the objects 
in the thin disk class. For the ADAF class objects, the X- 
ray spectra in the ROSAT and ASCA bands both steepen 
when TO increases. This contradicts the ADAF radiation 
model of the low accretion state with to ~ (0.01 — 0.08), 
where the calculated X-ray spectra should be harder and 
smoother because the optical depth goes up and causes a 
correspon ding increase in Compton y-param eter when m 
increases ( |Esin, McClintock fc Narayan 1997 ). However, it 
can be understood if wc adopt a model where the accretion 
disk consists of two zones: an outer standard thin disk ex- 
tending from a large radiu s down to rtr an d an inner ADAF 
from rtr down to r = 3 ( Narayan 19961 ). If rtr decreases 
with increasing m, and if rtr decreases below lO"'^'^, the 
ADAF X-ray emission spectra become dramatically softer 
because of the cooling effect when the radiation from the 



would be more important for thin disk and less important 
for an ADAF since the solid angle is larger than the former 
case. This would cause the X-ray spectra of the objects in 
thin disk class softer than those of the objects the ADAF 
class (Zdziarski, Lubinski, & Smith 1999). We can see in 
figures |l| and figures || that the spectral indices of the ob- 
jects in the thin disk class tends to be softer than those of 
the objects in the ADAF class, statistically. The accretion 
flows may have other forms, e.g. A DIOS (advection domi- 
nated inflow and outflow solution, Blandford & Bcgclman 



1999), which may be a possible alternative explanation for 



the relation. Detailed models should be explored to reveal 
more information on the two accretion classification. 
Recent observations show that 1H0419-577 is a "two- 



state" soft X-ray Seyfert 1 galaxy (Guainazzi et al. 1998). 
In its hard state, the entire 1.8-40keV spectrum can be 
well-described by a simple flat (F ^ 1.55) featureless 
power-law, while in its soft state, the soft X-ray spectrum 
is quite steep (T ^ 2.5). The 0.5-2.5keV flux changes by 
a factor of 6 from its soft state to its hard state. As- 
suming that Lion is proportional to the flux at 0.5-2.5keV, 
Lion/LBdd also changcs by a factor of 6. Based on its 
optical observation, the Eddington ratio of ionizing lumi- 
nosity is estimated to be less t han 0.08. So, 1H0419 -577 is 
an object in the ADAF class ( Guainazzi et al. 1998| ). The 
Log(Lion/LEdd) of this object is predicted to be ^ —1.72 
in its soft state by the fitted straight line for the ADAF 
class in section 3. Therefore, its Log(Lion/LEdd) is about 
-2.4 in the hard state. The inferred hard X-ray spectral 
index from figure || is consistent with the observations. 

We summarize our main result as the existence of two 
different accretion classes: the ADAF and the thin disk 
accretion class in Seyfert 1 galaxies and QSOs. The 
soft X-ray spectrum becomes steeper when the accretion 
rate increases in both classes. The detailed fitting of 
Seyfert 1 galaxy and QSO spectra with ADAF and/or thin 
disk models should be explored to reveal the underlying 
physics, as has been done for BHXBs and low-luminosity 
active nuclei. 

We thank an anonymous referee for useful comments 
which improve this paper. We thank Neta Bahcall, Scott 
Tremaine and Joseph Weingartner for a critical reading. 
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Fig. 1. — The ROSAT soft X-ray spectral indices versus the estimated Eddington ratio of ionizing luminosity. Filled symbols arc for radio- 
quiet objects; open symbols are for radio-loud objects; circles are for Seyfert 1 galaxies; triangles are for QSOs. The dashed line represents 
the critical Log (Lion/i'Edd) foi' the ADAF and the thin disk classes. The dot-dashed lines are the fitted lines for the two classes. 
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Fic;. 2. — The ASCA X-ray photon idices versus the estimated Eddington ratio of ionizing luminosity. The symbols are the same as in 
figure 1. The dashed line represents the relation of Fasca versus Log(Lion/i'Edd)- 



TABLE 1 

The black hole mass, the Eddington ratio of the ionizing luminosity and the X-Ray parameters. 





Name 


Redshift 


Mbh 


Log(- 




Rcf. 


Q^Rrm AT 


Rcf. 


Tasca 


Ref. 




WMq 


t'Edd 




3C 120 


0.033 


o 1 +2.0 


-0.57 


1 


1.15 + 0.33 


4 


1 


07+0. 06 
" ' —0.05 


13 




3C 390.3 


0.057 


oq 1+12.4 


-2.19 


1 
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13 
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a^ ■ 


Mrk 110 
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-0 
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Mrk 335 
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49 


1 


1.92 + 0.04 
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1 


•'"-0.13 


13 


> ■ 


Mrk 509 


0.035 




-0 


54 


1 


1.58 + 0.10 


4 


1 
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13 


o ; 


Mrk 590 


0.027 




-0 


89 


1 


1.41 + 0.06 


7 








Mrk 817 


0.033 


' -0.9 


-1 
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1.90 + 0.24 
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NGC 3227 
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0.009 
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10 
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NGC 4051 
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86 


1 
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9 
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13 
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NGC 4151 
NGC 5548 
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-1 9+0.8 
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49 
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1 
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4 
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i 


1 /I 

1 .40 
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13° 
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oo . 


NGC 7469 
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-1 
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1 
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4 


1 
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-0 
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Mrk 279 
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1 


QQ+0.08 

°'J-0.07 


13 




PG 0052+251 


0.155 


55.0 


-1 


11 


3 


1.27 + 0.12 


4 








o ■ 


PHL 909 


0.171 


372 


-1 


59 


3 












NAB 0205+024 


0.155 


3.55 


0.46 


3 


2.30 + 0.09 


4 


2 


-^"-0.09 


15 


a . 


PG 0923+201 


0.190 


123 


-1 


42 


3 










PKS 1004+130 


0.240 


126 


-1 


10 


3 












jv: 


PG 1012+008 


0.185 


11.7 


-0 


61 


3 


1.66+ • • • 


4 








HE 1029-140 


0.086 


200 


-1 


20 


3 


1.47 + 0.06 


10 










PG 1116+215 


0.177 


26.3 


-0 


43 


3 


1 68 ± 06 


11 










PG 1202+281 


0.165 


28.8 


-1 


34 


3 


1.23 + 0.11 


4 










3C 273 


0.158 


74.1 


-0 


30 


3 


1.07 + 0.08 


4 


1 


^q+0.01 

"J^-O.Ol 


14 




PKS 1302-102 


0.286 


56.2 


-0 


37 


3 


1.36 + 0.39 


11 








PG 1307+085 


0.155 


11.0 


-0 


45 


3 


1.64 + 0.16 


4 










PG 1309+355 


0.184 


14.1 


-0 


71 


3 


1.44 + 0.31 


10 










PG 1402+261 


0.164 


6.92 


-0 


27 


3 


1.89 + 0.20 


11 










PG 1444+407 


0.267 


14.5 


-0 


41 


3 


1.94 + 0.32 


12 










3C 323.1 


0.266 


148 


-1 


22 


3 


1.36 + 0.22 


11 










PKS 2135-147 


0.200 


123 


-0 


87 


3 


0.82 + 0.39 


4 










PKS 2349-014 


0.174 


117 


-0 


89 


3 













'^The averaged Fasca of multi-observations are adopted. 
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